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Анализируются условия, при которых в процессе обработки 
материалов резанием формируется хаотическая динамика. В 
ранее выполненных исследованиях в случае потери устойчиво- 
сти процесса резания в окрестности равновесной системы об- 
разуются предельные циклы или инвариантные торы. В отли- 
чие от этих исследований рассмотрен случай, когда свойства 
инструмента таковы, что за счет изгибных деформаций образу- 
ется нелинейная положительная обратная связь. Приводится 
математическая модель системы для этого случая. На основе 
цифрового моделирования с использованием пакета приклад- 
ной программы МАТЕГАВ исследуется влияние параметров 
динамической модели при условиях образования хаотической 
динамики. В результате выполненных исследований установ- 
лено, что по мере увеличения параметров, характеризующих 
образование положительной обратной связи, система претер- 
певает ряд бифуркаций удвоения периода в системе образова- 
ния странных аттракторов. Они расположены в окрестностях 
точек равновесия и имеют ограниченную область. Показано, 
что хаотичные колебания инструмента приводят к хаотичному 
формированию обрабатываемой поверхности, поэтому в при- 
кладной области необходимо выбирать параметры, при кото- 
рых хаотичная динамика не образуется. Несмотря на то, что 
рассмотренные примеры относятся к процессам обработки 
резанием, полученные результаты имеют общее значение для 
динамических систем, взаимодействующих с различными сре- 
дами, например, с трибологической средой. 


Ключевые слова: процесс резания материалов, динамическая 
система, инвариантные многообразия, хаотические аттракторы, 
бифуркации. 


Введение. Динамическая система процесса резания, базовая модель которой характеризует взаимодействие про- 
странственной динамической модели подсистемы инструмента с динамической связью, формируемой процессом об- 
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работки, может служить примером, иллюстрирующим различные эффекты нелинейной динамики. В известных иссле- 
дованиях в области динамики процессов обработки на металлорежущих станках главное внимание уделялось рас- 
смотрению устойчивости точки равновесия системы, а также автоколебаниям системы [1-18]. Во всех этих работах не 
принимался во внимание случай, когда за счет существенных изгибных деформационных смещений инструмента, 
вызывающих уменьшение переднего угла режущего инструмента, силы резания не уменьшаются, как в отмеченных 
выше работах, а возрастают. Тем самым формируется положительная обратная связь, способствующая самовозбужде- 
нию системы резания. Ниже будет показано, что в этом случае, как правило, в окрестностях равновесия образуются 
странные (хаотические) аттракторы. Приводимый в статье материал дополняет известные примеры образования хао- 
тических аттракторов [19—24]. 


Математическая модель. Как и в ранее опубликованных работах, будем рассматривать базовую динамическую мо- 
дель механической части системы в виде [11-14], (рис. 1): 
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ь Т 
мента в плоскости, нормальной к поверхности резания; Е\(Г) = {Е (Е), Р›(ГР)} — вектор сил резания. Упругие деформа- 
ционные смещения рассматриваются в подвижной части суппорта Ге. Угловая скорость вращения заготовки считает- 


ся постоянной. Она определяет постоянную составляющую скорости резания Г. Поэтому в установившемся состоя- 
нии при условии устойчивости положения равновесия величина подачи на оборот 5, = сопзЁ . Сила резания формиру- 


ы 1 
ется в виде суммы двух составляющих, действующих на переднюю грань инструмента К (и на его заднюю поверх- 


2 1 2 ь 
ность Е ) то есть Е=Е® + Е) (рис. 1). Однако в рассматриваемом случае влиянием сил, действующих на зад- 
нюю грань инструмента, можно пренебречь, то есть угол © есть величина большая, что характерно, например, для 


растачивания. Технологические режимы: глубина резания {, и величина подачи на оборот 5, при заданной геомет- 


рии инструмента однозначно определяют ширину срезаемого слоя Б и его толщину а, так как а=брзшф, а 


Ь=1,/заф (ф — главный угол инструмента в плане (рис. 1)). Деформационные смещения в направлении, нормаль- 


ном к плоскости (А — Х›) не принимаются во внимание, так как деформации в этом направлении практически не вы- 


зывают изменения площади срезаемого слоя э. Параметры т, Й и с можно определить по правилам, изложенным, 
например, в [11]. 

Главной особенностью рассматриваемого случая является изменение переднего угла режущего инструмента 
при деформационных смещениях инструмента в направлении Х,. В результате имеет место уменьшение значений 


переднего угла вплоть до отрицательных его значений. Поэтому, во-первых, по мере увеличения деформационных 
смещений в направлении Х, имеет место увеличение объема пластической деформации в зоне резания, и, как след- 
ствие, увеличение модуля силы резания [26]. Во-вторых, силы резания в направлении Х’ возрастают непропорцио- 


нально быстро [11]. Кроме этого, как и в ранее рассмотренных случаях [27, 28], имеет место запаздывание изменения 
сил по отношению к деформационным смещениям инструмента. Эти особенности процесса резания можно учесть, 
если рассмотреть силы резания в виде следующих зависимостей: 
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где Х\, Х, — деформационные смещения вершины инструмента в двух ортогональных направлениях в [мм]; 1» Х> 


2 2 
— угловые коэффициенты, удовлетворяющие условиям нормировки (х/)“ + (%>) =1; а, В — соответственно толщи- 


на и ширина срезаемого слоя, зависящие от технологических режимов и геометрии инструмента в [мм]; © — коэф- 


фициент, определяющий приращение силы к приращению объема деформируемого материала, имеющий смысл плот- 

1 2 1 2 
ности в [кг/мм’]; В — составляющая коэффициента жесткости процесса резания в [кг/мм|; р ), г т и ), т 
— постоянные времени, определяющие запаздывание изменения сил по отношению к деформационным смещениям 


инструмента относительно заготовки в [<]. Причем, (то и Т ов :(2) 





Рис. 1. Ориентация осей деформационных смещений и сил, действующих на переднюю грань режущего инструмента 


Таким образом, математическая модель динамической системы резания представляется в виде следующей 


системы: 
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которая учитывает все основные свойства динамической системы резания, отмеченные выше. 
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Свойства равновесия. Равновесие системы определяется исходя из условия: изменения деформационных смещений 


* * 
инструмента отсутствуют. Следовательно, точка равновесия (Х\,Х.) определяется из уравнения 


(са арЬ)Х! + с 1Хо = Х.рба + а(Х>)?; 


я : (4) 
(с12 +Х2РрБ) Лт + (2 -В)Х› = Хорва. 
Анализ (4) приводит к необходимости рассмотрения следующего уравнения: 
о(Х>) +С; Хо +Е, =0, (5) 


С. 2 (11 +Х1рБ) —[<11В + ХАРБВ + (12 ) + Х2рЬс,2]| А-В 
У До— Е — обобщенная жесткость системы; 


где С 
` (12 +Х2РБ) р 


161.2 — №26 ь 
= рра — суммарное значение силы. В динамической системе резания параметры А, В, О исум- 
(с12 +Х2РБ) 
марная сила Е; варьируются в достаточно широком диапазоне. Они определяются элементами матрицы жесткости 


подсистемы инструмента, свойствами процесса резания, технологическими режимами и геометрией инструмента, 
влияющей на угловые коэффициенты. 


Вначале рассмотрим случай, когда х.с! > =Хс! 1. Возможность выполнения этого условия зависит от геомет- 


рии инструмента, которая определяет угловые коэффициенты, и углов ориентации осей эллипсов жёсткости в подси- 
стеме инструмента, то есть от конструктивных свойств суппортной группы. Для этого случая наибольший интерес 


представляет выяснение влияния на свойства системы параметра В. Из (5) видно, что по мере увеличения В с учетом 


5 
(>) ., 
соотношения ус. =Х.с| существует точка бифуркации Во = с›2 -———_ с › начиная с которой имеет место ветв- 


(Х1) 


ление решений. Это бифуркация равновесия системы типа вилки. Ей соответствует С; =0. Как видно, при В\Ву си- 


стема имеет единственную точку равновесия (рис. 2, а). При В)Ву свойства системы, рассматриваемые в вариациях 
* 
относительно двух точек равновесия, например, + Х.1, являются симметричными, как это наблюдается в системе 
“ “ * 
Лоренца [19]. Не трудно показать, что для этой системы в потенциальной функции точкам + Х., соответствуют два 


<“ “ * 
минимума. При В‹Во потенциальная функция имеет единственный минимум, соответствующий точке Х> =0. 
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Рис. 2. Графическая интерпретация формирования точек равновесия системы: а — случай, когда со = Х2С1 1; 


Ь — случай, когда Со = Х2С1 1 
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В общем случае условие х1с!> = Ус! Не выполняется. Тогда точка бифуркации и точки равновесия системы, 


во-первых, становится зависящими от К; , во-вторых, свойства системы в вариациях относительно точки равновесия 
при неизменных параметрах также изменяются (рис. 2, Ь). В этом случае точке бифуркации равновесия соответствует 
условие С;)0. Отметим, что точкам равновесия Х> из (4) однозначно соответствуют Х!. 

Для анализа устойчивости точек равновесия необходимо рассмотреть линеаризованное в окрестности равно- 


весия уравнение (3), получаемое после замены Х/ (В) = Ху +х (6) и Х, (6 =Х. +х› (1) 
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а ао (6) 
4” @ 
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—. й,> — ХорЬТ! й > + Вт с1.2 +Х2РЬ с›,2 -В 


Для анализа устойчивости конкретной точки равновесия необходимо рассматривать корни характеристиче- 
ского полинома системы (6) 


тр” + (в -жРЬТ)р+сы +2аРЬ [пол +3 То (Хр + сл -З0(Х 5)" 
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А(р) = 


Анализ ЭВОЛЮЦИИ корней характеристического полинома показывает большое разнообразие возможных меха- 


низмов потери устойчивости равновесия системы. В частности, если варьируется суммарное значение силы Ёу ‚ ВЛИЯ- 


* < 
ющее на Х> ‚ То знак смещения точки равновесия не влияет на своиства системы в вариациях относительно точки 
авновесия. Мат ИЦЫ й и с; в (6) являются несимметричными, то есть п едставимы в виде сумм симметричных и 
У. > ) 


кососимметричных составляющих. Анализ конкретных систем показывает, что главным механизмом потери устойчи- 


вости равновесия связан с преобразованием симметричной составляющей матрицы й; из положительно определен- 
ной в отрицательно определенную. Кроме этого, по мере увеличения В возможно преобразование симметричной (по- 
тенциальной) составляющей матрицы с; из положительно в отрицательно определенную. При варьировании В и РЁ; 


имеют место множество бифуркаций многообразий, формируемых в окрестности равновесия. 
Анализ бифуркаций системы. Проанализируем бифуркации системы на конкретном примере. При этом будем учи- 
тывать, что зачастую изучение примеров позволяет выяснить общие свойства. Так как при анализе бифуркаций необ- 
ходимо рассматривать большие отклонения от точек равновесия, то для исследования необходимо рассматривать си- 
стему (3). Решение этой системы не может быть найдено в явном виде. Поэтому для ее интегрирования использовался 
метод Рунге-Кутта 4-го порядка. 

0,2.10—* 0 12 08] [800 900 


Рассмотрим две системы, имеющие параметры: т= 5 — . С = 
ь И 0 02.10? 0,8 2,0 900 1200 


Размерность элементов матриц соответственно равна: Кг. с? } мм], [кг-с/мм|, [кг/ мм]. Они отличаются угловыми 
коэффициентами ориентации силы резания. Для системы №1 — 9, =0,8, Хх, =0,6; для системы №2 — , =0,6, 
Хх = 0,8. Таким образом, для системы №1 по мере увеличения В ветвление равновесия симметрично, симметричны и 
свойства в вариациях относительно точек равновесия. Для системы №2 имеет место асимметрия. Параметры процесса 


2 4 
резания остаются неизменными и соответствуют: В =3,0мм, а=0,2мм, р=2000кг/ мм”, В=0,9.10`кг/ мм | 


т = 0,002с, 1 = 0,003с, 7 = 0,001е, 7? = 0,002с [12]. 

На рис. 3 приведены две диаграммы бифуркаций многообразий стационарного состояния в окрестности рав- 
новесия при варьировании Ве (0,2000)кг/ мм. Проанализируем их по мере увеличения параметра В. Вначале про- 
анализируем стационарные многообразия для системы № 1. По мере увеличения В единственная точка равновесия 
является асимптотически устойчивой (множество М1,В=[0-— 755] ), затем имеет место бифуркация точки равновесия, 


при этом верхняя и нижняя точки двух ветвей диаграммы равновесия являются асимптотически устойчивыми (множе- 
ство №2). Причем, устойчивость точек равновесия является апериодической. Затем (множество №3) наблюдается би- 
фуркация Андронова-Хопфа рождения орбитально асимптотически устойчивой пары предельных циклов. Причем эти 


эинэГэчонишеи и эинэодтоонише 


15 


БИр://уезыиК.Аоп$ва.га 


Вестник Донского государственного технического университета 2015, №2(81), 11-21 





предельные циклы являются симметричными в окрестностях верхнего и нижнего равновесий. Они располагаются в 
ограниченном пространстве, расположенном между двумя верхними и нижними ветвями. 
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Рис. 3. Бифуркационная диаграмма при варьировании В: а) 2; =0;) Е; =0 


При дальнейшем увеличении В (множество №4) наблюдается каскад бифуркаций удвоения периода. Извест- 
но, что бифуркация удвоения периода характеризует один из сценариев рождения хаотического аттрактора, рассмот- 
ренный М. Фейгенбаумом [25]. Хаотические аттракторы формируются на множестве №5. Причем они образуются в 
окрестностях верхней и нижней точек равновесия, существуя в ограниченном объеме фазового пространства. Одно- 
временно они характеризуют притягивающее многообразие. Далее, за счет увеличения интенсивности возбуждения, 
траектории становятся орбитальными относительно всех трех точек равновесия. При этом вначале формируется тра- 
ектория с каскадом удвоения периодов (множество №6). Затем после каскада бифуркаций удвоения периода вновь 
формируется хаотический аттрактор (множество №9). Однако, в отличие от ранее рассмотренного хаотического ат- 
трактора, траектории захватывают все три точки равновесия. В дальнейшем циклы преобразования стационарных тра- 
екторий периодически повторяются (множества №10, №11). Наконец, система становятся неустойчивыми в целом 
(множество №12). Однако в этом случае траектории, двигаясь относительно каждой точки равновесия, находятся в 


пространствах, которые не пересекаются друг с другом. Заметим, что, начиная с В=1500кг/ ммсвойства системы 


становятся чрезвычайно чувствительными к вариациям этого параметра. 


Некоторые наиболее типичные примеры проекций фазовых траекторий на плоскость Х, — АХ. / & приведены 


на рис. 4. Здесь необходимо обратить внимание на то, что хаотические аттракторы при ВеМ4 формируются исклю- 


чительно в окрестности двух точек равновесия, а при Ве М9 они напоминают аттрактор Лоренца и имеют свойства, 


близкие к этому аттрактору [23, 24]. 
Если Ру = 0 (система № 2), то бифуркационная диаграмма за счет асимметрии верхней и нижней ветвей явля- 


ется более сложной. По мере увеличения |} система вначале также имеет единственную точку равновесия, которая 


является асимптотически устойчивой (множество М1). Затем, в окрестности единственной точки равновесия форми- 
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руется устойчивый предельный цикл (множество №2), бифуркация удвоения периодов (множество №3), затем хаоти- 
ческий аттрактор (множество №4). Процесс удвоения периодов (множества №5 и №7) и формирование хаотических 


аттракторов (множества №6 и №8) повторяется вплоть до точки, в которой наблюдается ветвление равновесия 
Во =1386кг/ мм. 
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Рис. 4. Некоторые примеры проекций стационарных траекторий на фазовую плоскость Х›, АХ. / ЧЁ по мере увеличения 


параметра В для Е; =0:а)№2;5) №3; ©) №4; 4) М5 ;е №6; Л №9 


Затем свойства системы в окрестностях верхней и нижней ветви равновесия существенно меняются. Вначале 
в окрестности верхней ветви формируется странный аттрактор, а нижней — асимптотически устойчивая точка равно- 
весия (множество №9). Затем в окрестности верхней ветви (множество №10) в системе вновь происходит бифуркация 
удвоения периода циклов, а на нижней — устойчивая точка равновесия. Наконец, в области нижней ветви наблюдает- 
ся бифуркация Андронова-Хопфа (множество №11), в области верхней ветви продолжается процесс удвоения перио- 
да. При дальнейшем увеличении В в области верхней ветви формируется хаотический аттрактор, а в области нижней 
ветви — сохраняются автоколебания (множество №12). Как предельный цикл, так и хаотический аттрактор характе- 
ризуют притягивающие стационарные многообразия, существующие в ограниченном фазовом пространстве. Наконец, 
в области верхней ветви система теряет устойчивость, и траектории уходят от центральной ветви, в области нижней 


ветви в ограниченном пространстве сохраняются автоколебания (множество №13). При дальнейшем увеличении В 


система в области верхней точки равновесия остается неустойчивой, а в области нижней ветви в ограниченно про- 
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странстве вначале образуется каскад бифуркаций удвоения периода (множество №14 ), а затем образуется хаотический 
аттрактор, занимающий ограниченную область в окрестности нижней точки равновесия (множество №15). За преде- 
лами этой ограниченной области система движется по траекториям, отходящим от всех трех точек равновесия. Нако- 
нец, (множество №16 ) система становится неустойчивой в целом. 


Принципиальным отличием каскада бифуркаций при РА; *0 от систем Е; =0 является различное поведение 


ее в областях верхней и нижней ветвей траектории смещения точек равновесия (рис. 5). При этом не образуются ста- 
ционарные траектории, включающие сразу три точки равновесия, как это наблюдается в системе № 1. Обратим вни- 
мание на возможность существования в неустойчивой в целом системе некоторых локальных областей (рис. 5, с, @), в 
которых система является устойчивой в малом (рис. 5, с). Существуют также случаи, когда формируются хаотические 
аттракторы в ограниченном объеме (рис. 5, 4). Цифровое моделирование показывает, что на формируемые многообра- 
зия, кроме В и Е; ‚ оказывают влияние практически все параметры системы (3). Если проранжировать параметры по 
чувствительности к ним изменения многообразий, то необходимо отметить существенное влияние постоянных време- 
ни в (3), которые моделируют запаздывание сил по отношению к деформационным смещениям. Именно через меха- 
низм влияния запаздывающих аргументов объясняется формирование многих стационарных многообразий, рассмот- 


ренных выше. 





Х,,мм 107 


с) 4) 





Рис. 5. Пример преобразования траекторий вершины инструмента в плоскости Х| — Х› 
по мере увеличения параметра В для Р; =0:а) №11; 5) “12; ©) 514; а) “15; е) “16 
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Выводы. Динамическая система резания может служить примером сложных преобразований стационарных многооб- 
разий, формируемых в окрестности точек равновесия. В зависимости от параметров динамической характеристики 
процесса резания пространство состояния динамической системы резания структурируется, и при этом свойства си- 
стемы становятся чувствительными к малым вариациям параметров. Характерными для системы являются бифурка- 
ции Андронова-Хопфа, удвоения периода колебаний и образование хаотических аттракторов. В рассматриваемой си- 
стеме не обнаружено формирование многообразий типа инвариантного тора. Необходимо подчеркнуть, что все эти 
преобразования имеют практическое подтверждение, экспериментально полученное нами и другими исследователями 
[11, 29]. Если принять во внимание, что система резания является возмущенной, то в ней в реальных условиях образу- 
ется сложный, трудно предсказуемый режим колебаний. Важно подчеркнуть, что эти колебания оказывают влияние 
на параметры качества изготовления деталей. 
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